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We investigate the DD production in e + e~ annihilations near threshold in an effective Lagrangian 
approach. This shows that the lineshape of the cross section near threshold is sensitive to the 
contributions from 0', though it is below the DD threshold. The recent experimental data from the 
BES and Belle collaborations allow us to determine the ip'DD coupling constant which appears to 
be consistent with other theoretical studies. As a consequence, the 0'-0(377O) mixing parameter 
£NJ ■ can be extracted around the 0(3770) mass region. Resonance parameters for 0(3770), X(3900), 

0(4040), and 0(4160) are also investigated. The X(3900) appears as an enhancement at around 3.9 
GeV in the Belle data. In addition to treating it as a resonance, we also study the mechanism that 
the enhancement is produced by the DD* + c.c. open channel effects. Our result shows that such a 
possibility cannot be eliminated. 

PACS numbers: 13.66.Bc, 12.38.Lg, 14.40.Gx 
^ , I. INTRODUCTION 

(f) , As the first charmonium state above the DD threshold, the production of 0(3770) in e + e~ — > DD 

serves as a peculiar probe for exploring the QCD dynamics in the interplay of the perturbative and the 
non-perturbative regime. During the past few years experimental measurements were performed in the 
0(3770) region and some interesting observations were exposed. First, BES-II reported a branching ratio 
for 0(3770) — > non-DD up to 15% [HQ, while CLEO-c found a much smaller non-DD branching ratio Q. 
Note that so far only one exclusive channel -0(3770) — > (jyrj has been observed in 0(3770) non-charmoniun 
strong decays. One major concern is how the light hadrons are produced in cc annihilation. In other 
words, whether this process is dominated by pQCD or there are signs for non-pQCD contributions are 
questions for both experiment and theory. 

Further interests in this issue were raised by a recent non-relativistic QCD (NRQCD) calculation to 
the next-to-leading order (NLO) for the ■0(3770) non-DD decays Q, where the authors found significant 
QCD corrections from NLO. Implications of such a result would be that non-pQCD mechanisms may 
start to play a role. Quantitative studies of non-pQCD mechanisms are presented in Ref. @, where 
the authors show that the long-distance interactions due to intermediate D meson loops are essential for 
understanding the non-DD decays of the -0(3770). This also suggests that a dynamic understanding of 
the correlation between the 0(3770) non-DD decay and the so-called Okubo-Zweig-Iizuka rule evading [7[ 
mechanism is required. 

■0(3770) production in e + e~ annihilation is also useful for further studying the properties of the 0(3770). 
As shown by BES measurement an obvious deviation of the resonance excitations from Breit-Wigner 

is observed in e + e _ — > DD. This raises the question about the sources of causing such a deviation, and 
the role played by background processes [10]. One possibility is that some new structures in the energy 
region between 2.70 and 3.87 GeV may cause such a lineshape anomaly js|. Above the 0(3770) mass, 
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the data from Belle Collaboration suggest an enhancement around 3.9 GeV [TT|, which could be a signal 
for resonance. Note that 3.9 GeV is at the open channel threshold for DD* + c.c. Thus, it would be 
interesting to investigate the DD* + c.c. open channel effects and compare them with the Breit-Wigner 
solution in the numerical fits. By clarifying these issues in both experiment and theory we expect that 
the 0(3770) resonance parameters can be determined. Furthermore, the recent controversial results from 
BES-II [ISI and CLEO-c [|| on the "0(3770) non-DD decays can be disentangled [lEI. 

In this work we will study the lineshape of the e + e~ — > DD cross section in an effective Lagrangian 
approach. We will show that the cross section experiences important interferences from the ip' which will 
account for the ■0(3770) lineshape deviation from the Breit-Wigner form. A recent analysis of this issue 
was also done by Ref. [l2|]. In our work, we will show that the near threshold cross sections not only 
provide evidence for the ip'-ip (37 '70) interferences, but also provide a peculiar constraint on the dynamics 
for ip' couplings to DD. By determining the g^ijjD coupling, we can then examine the ipi/2SY]b(lD) 
mixing for the ip' and ^(3770), and extract the mixing angle at the mass of the V(3770) Iall3j. We 
will study the energy-dependence of the mixing parameter. It is essential to keep unitarity, and further 
dynamical information about the state evolutions could be gained. 

We will also investigate the bump around 3.9 GeV observed by Belle [111 ]. Given the success of the 
potential quark model (see Ref. [l4| for a recent review), a vector charmonium with J PC = 1 at 
3.9 GeV will cause great concern about the non-relativistic cc phenomenology. We will show that this 
enhancement may be caused by the DD* + c.c. open channel effects. 

This work is organized as follows: In Sec. II, we formulate the transition amplitudes for e + e~ — > DD 
from different sources, which include charmonium resonance excitations, electromagnetic (EM) back- 
ground from the vector meson dominance (VMD) model and open DD* + c.c. effects via intermediate 
meson loops. In Sec. Ill, we present our numerical results along with experimental observables. A 
summary will be given in Sec. IV. 

II. THE MODEL 

The ingredients considered in this work include charmonium resonance excitations and EM background 
(Fig. [T]), and open-charm effects via intermediate meson loops (Fig. [5]). The charmonium resonance 
excitations are constructed in the VMD model [IB - [l7l |. where the EM field will be decomposed into 
vector meson fields with both isospin and 1 components. By taking away explicitly the nearby resonance 
contributions from e.g. "0(3770) and 0', the EM background contributions arising from the continuum 
part can thus be parameterized by an effective coupling and minimized in the numerical fit. 

We also examine the DD final state interactions in e + e~ — s- 7* — s- V' -> (DD, DD* + c.c, D*D*) — > 
DD due to intermediate DD, DD* +c.c, and D*D* meson loops as shown in Fig.[5J These contributions 
do not double-count the ip' excitation in Fig. [1] since a Breit-Wigner is explicitly introduced there. 
The effects of such a loop contribution may give rise to both energy dependence of the tp' total width 
and a relative phase to the QED diagram. However, we note in advance that these contributions are 
negligibly small. It also should be clarified that there is no need to consider the intermediate meson loop 
contributions for 0(3770) since the effective coupling 3^,(3770)01) i s extracted from experimental data and 
should have included the loop effects. 

To study the DD* + c.c. open channel effects, we apply the experimental data for e + e~ — > DD* + c.c. 
from Belle [18]. This is equivalent to isolating out the DD* + c.c. open channel contributions from all 
possible sources, and can be compared with the Breit-Wigner fit. 

In order to evaluate the diagrams of Figs. [T] and [2j the following effective Lagrangians are needed: 

£vdd = g VDD {Dd^D-d^DD}V», 
Cvdd* = ~ig VD D^a0^d a V p d»D*"D + h.c., 
C vd *d* = ^g V D*D*t a i3^d a V* p d^D* u T + h.c., 
Cvdd' = g D , r5 {Dd fl V-dt l DV}V*i*, 

Cvd'D* = g VD * D ^[(D*»dlD: - D*%Dt) - D*»dtD; + D^Km , (1) 
where V and V are the pseudoscalar and vector mesons, and e Q/ 3 M „ is the antisymmetric tensor. 
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A. Charmonium excitations in VMD and EM background contribution 

According to VMD [HEH, the EM current can be decomposed into two parts. One is the hadronic 
part containing a complete sum over all isospin-0 and 1 vector meson fields, while the other is the so- 
called "bare photon" field. An empirical role played by the EM "bare photon" field is to assure the proper 
normalization of the physical photon field. Its contribution is generally small. Therefore, this part is 
minimized in the VMD model. 

In reality, it is not possible to include all hadronic vector meson amplitudes. A commonly adopted 
method is to include the vector meson resonances in the vicinity of the considered kinematics, and then 
treat the unknown part as an EM background which now includes the off-shell contributions from those 
faraway resonances and the "bare photon" amplitude. It would also be our strategy here to apply the 
VMD model to e + e~ — > DD near threshold. Figure. QJa) and (b) demonstrate the hadronic contributions 
and the EM background, respectively. 





(b) 

FIG. 1: Feynman diagrams based on the VMD model for e + e~ — > DD near threshold. Diagram (a) is for hadronic 
vector meson contributions, while (b) represents the EM background. 

Setting m e ~ 0, the coupling of vector mesons with the virtual photon can be extracted from experi- 
mental data by the VMD model QJ, [13 : 



e 

fib 



3F 



ip— fe+ e 



1/2 



(2) 



where \p e \ is the electron three- momentum in the vector meson rest frame, r^_,. e + e - is the partial decay 
width, and a e = 1/137 is the fine-structure constant. 

The amplitude of Fig. HJa) due to ip resonance excitations can be written as: 



T a = e v{k 2 )"/fj.u(ki)- 



1 mi 



1 



s s- mi + im^T,, 



(3) 



where is the total width of the charmonium resonance, and u{k\) and v{k,2) are the Dirac spinors for 
the electron and positron, respectively. The four-vectors pi and p 2 are the momenta for the final state 
D and D meson. 

We parameterize the EM background contributions to e + e _ —> DD as follows: 



1 



T b = e v(k 2 )7fj.u(ki)-g c (s)(pi -p 2 ) M , 



(4) 



where g c is an effective coupling of the EM background contribution to DD. As pointed earlier, this 
amplitude contains the "bare photon" contribution and contributions from other vectors which are not 
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explicitly included. In this sense, the value for g c will depend on how many vector charmonium states are 
included in the fitting. The inclusion of more higher states near threshold will minimize this amplitude. 
Combining Eqs.([3|) and (E]), one obtains the total amplitude for e + e~ — > DD 

_ e 2 v{k 2 ) lll u{k 1 ){p 1 - p 2 y r . . rn\, g^, Dt) \ - g^. DB ., 

— \9c\ s ) > j ~2 — r~ — f — h Z,i 2 — — f — e J' 

s /^,/ s - m^, + im^'L ^ 7V< s - + im^l ^ 

where the ?// amplitude is explicitly included, and a phase factor e l ^ i is added to other charmonium 
resonance amplitudes above the DD threshold. Then, the cross section for e + e~ — > DD can be written 
as 

/ + - in m 87rQ e \P\ 3 I / \ j m l' , V"^ g^DD i^|2 (R s 

a{C C -+ DD) = —^ l9cis) + J^ S - m%, + i^IV + ^ - < + 6 ' ' ( ) 

where is the D-meson three- vector momentum in the overall cm. frame. 



B. Intermediate meson loop contributions and DD* + c.c. open channel effects 

The intermediate meson loops via the DD, DD*+c.c. and D*D* rescatterings in Fig.[5]can be evaluated 
as follows [EQIi: 

T L = -te v(k 2 h ll u{k 1 )-—^ 2— — / -— > . (7) 

s J j, s - ml + J (27Tj 4 ^ aia 2 a 3 

r ^ polarization 

Taking DD(p) as a example, the vertex functions for the DD{p) loop are 

T\ = igi(h - h) ■ £</, 

T 2 = ig 2 (h +pi) -£ P (8) 
T 3 = igz(h+V2) ■ £* p 

where g\, g 2 , and g% are the coupling constants at the meson interaction vertices (see Fig. [5]). The four- 
vector momentum, l\, h, and ^3 are for the intermediate mesons, respectively, while ax = if — ra\,a 2 — 
l\ — m|, and 03 = if — m\ are the denominators of the propagators of intermediate mesons. 

Divergences are inevitable in the loop integrals. Since the effective Lagrangian approach is not a 
renormalizable theory, a common way to kill the divergences is to introduce a form factor as a cut-off for 
those unphysical contributions from the ultraviolet regime. Such a prescription will also compensate the 
internal particle off-shell effects. The form factor JF(p 2 ) is usually parameterized as 

<*> 

where n = 0, 1, 2 correspond to different treatments of the loop integrals. If the on-shell approximation 
is applicable [20h22| ]. the loop integrals will not suffer from the divergences, and the form factor will take 
care of the off-shell effects. In the present work, we consider the dipole form factor, i.e. n = 2, in the full 
loop integral. The cut-off energy A is usually parameterized as 

A = m ex + oAqcd, (10) 

where Aqcd — 220 MeV, a is a tunable parameter, and m ex is the mass of the exchanged meson. There 
are different expressions for the form factor adopted in the literature. However, note that the cut-off is 
always tunable. We would generally need experimental data to determine a proper form factor parameter. 
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The charmed meson couplings to the light mesons are obtained in the chiral and heavy quark limits [2c 



2 i gD'Dir 

9D*Dtt = -r- 9V m Dm D ', 9D'D*tt = — = — ■ 
Jtt M D 

9D'D P = V2\g pi gDDp = 9d*d p Md, 



(11) 



with f n = 132 MeV, and Md = ymomH* is a mass scale. The parameter g p respects the relation 
9p = rrip/U H3. We take A = 0.56 GeV" 1 and g = 0.59 [H . 





FIG. 2: Feynman diagrams for intermediate meson loops as corrections for the DD coupling to a charmonium 
state. 




FIG. 3: Feynman diagram for the DD* + c.c. open channel effects 

The DD* + c.c. open channel effects (see Fig. [3]) can be calculated in a similar way. 
data from Belle for e + e~ — > DD* + c.c. [l8[ will allow us to extract the ^*DD* coupl 



ilar way. The experimental 
DD* + c.c. will allow us to extract the ^*DD* coupling g^ D o* (s) (l9l j: 

(12) 



a(e + e 



Air lr)1 3 

DD*+c.c.) = Y ^al\g rDl3 *(s)\ 2 



where \p\ is the D-meson three-vector momentum in the overall cm. frame. 

As briefly mentioned earlier, explicitly adopting the experimental data for e + e~ — > DD* + c.c. means 
that all the resonances which can couple to DD* + c.c. are included in the effective coupling form factor 
9-y*nn» ( s )- This will lead to double-counting with the treatment of explicitly including V>(3770) and other 
resonances. Because of this ambiguity we will only calculate the exclusive open channel contributions in 
the following section, and discuss its behavior around the open channel kinematics. 



C. ^'-^(3770) mixing 

The nonvanishing g^DQ extracted from the cross section measurement will allow us to investigate the 
ip'-ip (3770) mixing. As a dynamic consequence of the intermediate DD loop (and DD* + c.c. etc), such 
a higher order effect can be quantified in e + e~ — > DD. 



G 



For two-state tp' — -0(3770) mixing, the covariant propagator can be expressed as [26l |2 

Q _ \ [ D <1>' ^'^(3770) | ^ 13 \ 

- D ^'- D 0(377O) - 1^^(3770) I 2 V ^''^(3770) A/,(3770) /' 

where D^i and -0^,(3770) are the denominators for the propagators of ip' and ■0(3770): 

D i = «-m? + *Vir < («), (14) 

with the energy- dependent width 

r ' i(s) = e^^+u- fo+(s)|3 + 4^°5oIpdo( s )| 3 + ••■] + r^- DD , (is) 



dominated by the DD channel near the 0(3770) mass region. |pj(s)| = ys — imf/2 is the three-vector 
momentum carried by the intermediate D meson at an energy s/s. 
-0(3770) i s the mixing term via the DD meson loop: 

A/>'V(3770) = 6^ ^ l9ip'D+D-9ip{3770)D+D- \PD+{s)\ 3 + g^> D° D°9^(3770)D° D° \PD° {s)\ 3 H ] , (16) 



from which we can define the mixing parameter |£| 

161 



A/>'i/>(3770) 



(17) 



Several points should be made: (i) In the energy region near the 0(3770) mass (below the DD* + c.c. 

open channel), the DD loop is the dominant contribution to the mixing matrix element. We restrict 

the discussion in this region and neglect the contributions from the DD* + c.c. and other high threshold 

loops. But we mention that contributions from the DD* + c.c. loop are rather small when y/s is above 

the threshold, (ii) In principle, the vertex couplings g^ipD and 3^,(377o)d£> should be s-independent in 

Eq. ©. In our case, the g^i^jj is extracted by fitting the cross section data, and no obvious s-dependence 

is needed within the energy region near threshold. We thus adopt a constant coupling here. Nevertheless, 

it shows that the s-dependence of the total width r,(s) is also weak. This justifies that the mixing 

parameter extracted above is at an accuracy of leading order, (iii) In Eq. (fTS! the non-Dl) channel 

contributions to the total width are rather small due to much weaker couplings. As a leading order 

estimate we apply the PDG values to fix these two widths and do not consider their s-dependence, i.e. 
r non-DB = 31? KeV and r non-£>fl = 415 MeV ^ 



III. NUMERICAL RESULTS 



Proceeding to the numerical calculation, our fitting strategy is to minimize the background contribu- 
tions given that the necessary resonance amplitudes are included. We thus perform two separate fits of 
the data. One is restricted to the near threshold region (Fit-I) and the other is extended to yfs ~ 4.3 
GeV covered by the Belle data [ll| (Fit-II). This separate treatment is based on the different experimen- 
tal situations between BES and Belle. BES has a relatively detailed scan over the energy region from 
threshold to the upper end of the 003770) mass, while the energy bins in Belle data are much larger and 
the cross sections have larger uncertainties. Moreover, the Belle data are extracted from the exclusive 
initial state radiation (ISR) production of DD at cm. energy 10.58 GeV, and are corrected by the ISR. 
The advantage of the Belle data is that they cover nearly all the energies from threshold up to ~ 5 GeV. 



A. Fit-I 



Quantum mechanically, the most important contributions to the e + e — > DD cross sections near 
threshold would come from the nearby vector charmonia, i.e. 0(3770) and 0'. Although 0' is below the 
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TABLE I: Fitting parameters for e+e" -> D°D° and e+e" -> U+D" near threshold. The experimental data are 
from BES 0. 





e+e" -> 


e+e" D+£>~ 


9il>'DD 


9.05 ± 2.34 


7.72 ± 1.02 


9ip(3770)DD 


13.58 ±1.07 


10.71 ± 1.75 


"1^(3770 (MeV) 


3774.0 ± 1.3 


3774.0 ± 1.6 


1^(3770) (MeV) 


28.4 ±2.9 


29.6 ±2.9 




228.1° ± 18.6° 


190.6° ± 21.1° 


x7d.o.f 


12.28/9 


13.32/9 



DD threshold, its mass is not located far away. Thus, its off-shell contribution may still be sizable. Since 
the -0' has a well-defined mass position and total width, we then fix these two quantities in the analysis, 
but leave its coupling to DD to be fitted by data. 

Other fitting parameters include the mass and total width of -0(3770), the -0(3770)1)1) coupling, and 
the relative phase angle between the -0(3770) and -0' amplitude. By fitting these parameters to the 
BES data Q , we actually minimize the contributions from the other resonances and "bare photon" . In 
Table. HI the fitting results for e + e~ — > D°D° and e + e~ — > D + D~ are listed separately. It shows that 
the values of x 2 are compatible in these two channels. Interestingly, the extracted couplings g^ijjjj and 
9ip(3770)DD appear to have significant differences in the charged and neutral channel. As a comparison, if 
we extract g^ 7 o)DD b Y the PDG data [28], we obtain 5</,(3770)d° d« = 12 - 43 and g^ 770 ) D +D- = 12.89, 
which are apparently different from the fitted values. In particular, the BES data suggest a relatively 
larger g^fsY70)D°D coupling than g^(377o)D+D- ■ We note that in both the charged and neutral channel, 
the fitted mass and total width of -0(3770) are consistent with each other. It is also interesting to note 
that the extracted values for g^uQ are within the range from other approaches [29T432T ] . 

In Fig. HI the fitting results for both channels are presented. The short-dashed curves are exclusive 
■0(3770) cross sections while the long-dashed curves are exclusive 0'. The solid curves are given by the 
interfering amplitudes between -0(3770) and 0' with the relative phase angle = 228.1° and 190.6° for the 
neutral and charged channel, respectively. Apart from the coupling differences, the differences between 
the phase angles also significantly contribute to the interferences. 

As an estimate of higher resonance effects, we take the PDG average of the total width and the upper 
limit of the branching ratio Br(V>(4040) ->■ DD) < 0.2% 0] to extract 5^ (4 o4o)d5 ~ 0-34. The -0(4040) 
decay constant can also be extracted from the experimental data by Eq. ([2]), which gives e//^,( 4 o4o) = 
9.35 x 10 -3 . In Fig. 0] the exclusive cross sections from -0(4040) are presented as the dotted curves, 
which turn out to be negligible near threshold. This justifies the treatment that the cross sections from 
threshold to the upper end of the ^(3770) mass are dominated by ^(3770) and -0' interferences, and 
the background contributions (here, they are referred to the "bare photon" and higher resonances) are 
minimized. 

We also calculate the energy-dependence of the cross section ratio between e + e~ — > D + D~ and e+e~ — > 
D a D a in Fig. [5] to compare with the BES data. This shows an overall agreement, except that we should 
note that the ratio may start to deviate from reality at higher energies since other resonances and 
mechanisms would start to have stronger interferences. Since the cross sections appear to be sensitive to 
tiny discrepancies between the charged and neutral channel, a more precise measurement of this quantity 
would be useful for understanding the underlying mechanisms. 

B. Fit-II 

In this subsection, we fit the experimental data from Belle fTTI j by including higher resonances. Apart 
from V>(4040) and -0(4160) from PDG (28|, we also include a new state X(3900) in the fit. As shown by 
the Belle data for e + e~ —> D°D° [ll|, an enhancement is observed at about 3.9 GeV which could be 
signals for a charmonium state. Although we fit the data using a Breit-Wigner here, we shall argue later 
that this enhancement may be due to the DD* + c.c. open channel effects. 
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FIG. 4: The cross sections for e + e~ — > D°D° (left panel) and e + e~ — > D + D~ (right panel) fitted by the tp' and 
V>(3770) interferences (solid line). The short-dashed line stands for the exclusive 0(3770) cross section, while the 
dashed line denotes the one for ip' . The dotted line represents the exclusive contribution from -0(4040). The data 
are from BES 

1.00-r 

0.75- 



0.50- 
Q£ ■ 
0.25- 

0.00- 

3.74 3.75 3.76 3.77 3.78 3.79 3.80 

Vs 

FIG. 5: The cross section ratio between e + e _ — > D + D~ and e + e~ — > D°D° versus cm. energy yfs. The data 
are from BES @|. 




The strategy is similar to Fit-I. Namely, we fix the mass and total width of the ip' and then leave 
9t(,>D"D to be fitted by data. The other fitted quantities include the resonance parameters for 0(3770), 
X(3900), 0(4040), -0(4160), and their relative phase angles to the 0' amplitude. In Table. HEl the fitting 
results are listed and can be compared with Table. HI 

This shows that the coupling g^jjopo appears to be stable in this fit, while parameters for ?/)(3770) 
exhibit strong sensitivities to the data. In particular, the total width becomes rather small, even smaller 
than the PDG average [28j . The phase angle also changes, and the overall effects are that the cross sections 
at the -0(3770) mass are increased. Such a dramatic change is likely due to the significant discrepancies 
between the BES and Belle data. Again, a detailed scan over the threshold region is strongly required. 
The fitted total widths for X(3900) and 0>(4O4O) turn out to have large uncertainties. 

We plot the fitting results in Fig. [51 which appear to be in good agreement with the Belle data. In 
comparison with Fit-I, the importance of -0' is consistently highlighted. Exclusive cross sections from 
other resonances are also presented in Fig. [6] We note that the smooth-out behavior at the high energy 
end is favored by the minimization, and it leads to a negligibly small contributions from 0(4160). Note 
that the data still have rather large uncertainties. The contributions from -0(4160) should be restudie 
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TABLE II: Fitting parameters from Fit-II. The results are obtained by fitting the Belle data [Tll | for e + e — > 
D°D°. We separately list the other two fitted quantities here: g^,i£,or,a = 9.05±0.37, g^(z770)D°D a ~ 13.13±0.65. 
The reduced % 2 is xV d -°- f = 17.45/19. 





M (GeV) 


T (MeV) 


g V D D /fv(xl0-*) 


4> 


X(3900) 


3.894 ± 0.011 


89.8 ± 12.6 


6.76 ± 0.89 


104.36° ± 7.90° 


V(3770) 


3.7724 ± 0.002 


25.4 ± 1.4 


23.2 ± 1.1 


198.85° ±4.19° 


^(4040) 


4.0812 ± 0.008 


96.2 ± 11.4 


3.48 ± 0.34 


101.16° ± 11.742° 



with high-quality data, and we will not discuss them much here. 




3.7 3.8 3.9 4.0 4.1 4.2 4.3 

Vs 

FIG. 6: The Belle cross section for e + e~ — > D°D° 11] fitted by Fit-II. The solid line represents the overall 
results, the long-dashed line denotes the exclusive contribution from ip' , and the long-dotted, short-dashed and 
short-dotted lines denote the exclusive cross sections from f/>(3770), X(3900) and V'(4040), respectively. 



C. The DD* + c.c. open channel effects 

The enhancement at 3.9 GeV from the Belle data, if is due to resonance excitation, may cause some 
confusion about the potential model predictions (l^ | for the charmonium spectrum since one would not 
expect a J = state in this mass region. Taking into account the success of the potential quark 
model, an alternative explanation for this structure would be due to the DD* + c.c. open channel effects 
as illustrated by Fig. [21 

With the help of Eq. (|T2|) and experimental data for e + e~ — > D + D*~ +c.c. from Belle [18], the effective 
coupling g 7 *£)5*(s) can be extracted. The corresponding datum points are plotted in Fig. [7] in terms of 
x (= s - (m D + tod-) 2 ). 

We employ two methods to extract the effective coupling form factor: 

(i) Form factor one (FF-I) 

In the FF-I scheme, we fit the data in Fig. [7] with an exponential function: 

9~f*DD*( s ) =3i ex P[-( s - ( m D +m D ,) 2 )/t 1 ] +g , (18) 

where x = corresponds to the DD* +c.c. threshold, and g\, t\, and go are fitting parameters. The above 
function can perfectly describe the energy dependence of the effective coupling. By simply extrapolating 
the exponential to the DD threshold (dashed line in Fig. [Jj, we then apply this form factor to calculate 
the open channel cross sections. 

(ii) Form factor two (FF-II) 
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TABLE III: Parameters fitted with FF-I and FF-II for the 7* DD* effective coupling form factor. The data are 
from Belle [3. 



FF-I 


FF-II 


9fDfi« = ffiexp[-(s - (mr> +m D -*) 2 )/t 1 ] + g 


9~<»DD* \ s -m\+im X T x 1 &1 


9i 


8.86 ±0.59 


fee. 


3.08 ±0.31 


ti 


1.28 ±0.06 


m x (GeV) 


3.943 ±0.014 


go 


0.43 ± 0.02 


Fx (MeV) 


119.0 ± 10.0 




61 


0.016 ±0.045 


xVd.o.f 


68.86/53 


xVd.o.f 


47.67/52 



In the FF-II scheme, we fit the data in Fig. [7] with a single resonance: 



9j'DD*{s) 



bo 



im x T x 



(19) 



with a background term b\. The parameter bo can be regarded as the product of the j*X coupling and 
the XDD* coupling. This parametrization agrees with the data at higher energies, but drops at the 
threshold region. 

The above effective couplings should be taken with caution. In principle, the parametrization of the 
FF-I scheme contains all the contributions from the vector resonances, while the FF-II scheme contains 
only one resonance plus a background. Comparing these two form factors with each other would suggest 
the dominance of one single resonance around 3.94 GeV, which seems to support the need for A" (3900). 
However, it should be recognized that the Breit-Wigner structure fitted by the form factor is largely due 
to the collective contributions from nearby resonances such as ^(4040). The parametrization somehow 
tilts the resonance signal due to the unphysical background term. In fact, both Belle [18| and BABAR [33j 
observe the dominant ■0(4040) in e + e~ — > DD* + c.c, but there are no signs for the A"(3900) there. This 
could, in fact, support the fact that the enhancement at 3.9 GeV in e + e~ — ¥ DD is due to non-resonant 
mechanisms. For our purpose of extracting the effective j*DD*(s) form factor, these parametrizations 
should be acceptable. 



> 





FIG. 7: The ^-dependence of g^-oD" extracted from the cross sections for e e — > D D* + c.c. [18| with 
x = s — (mo ± m_o*) 2 . The dashed and dotted line are given by the FF-I and FF-II schemes, respectively. The 
vertical line at x = labels the DD* threshold. 



With the form factor coupling g 1 *oD-" ( s )i we can then evaluate the DD* + c.c. open channel effects 
in e + e~ — > DD. We plot the cross sections from the DD* open channel in Fig. [5] As a comparison, we 
also include the Breit-Wigner form factor (solid line) for AT(3900). Taking the cut-off parameter a = 2.5, 
the dashed and dotted line are given by those two form factors, FF-I and FF-II, respectively. This shows 
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3.7 3.8 3.9 4.0 4.1 4.2 4.3 



Vs 

FIG. 8: The comparison of cross sections from different mechanisms around 3.9 GeV. The solid line is given by 
a Breit-Wigner resonance X(3900) in Tab. [TTJ while the dashed and dotted line are given by the FF-I and FF-II 
schemes, respectively. The vertical line labels the DD* threshold, a — 2.5 has been taken. 

that FF-I can lead to some enhancement to the cross sections below ~ 3.9 GeV, while the cross sections 
from these two form factors are similar to each other above 3.9 GeV. This shows that the behavior of 
the solid line is similar to FF-II except that the fitted mass and width are slightly different. The reason 
again should be that the fitted form factor FF-II contain all contributions from several nearby states. 
Thus, the fitted Breit-Wigner parameters contain unphysical information. 

It is interesting to see that the cross section peaks from these different treatments are located in a 
similar place. In this sense, we argue that the present data cannot eliminate the possibility that the cross 
section enhancement at around 3.9 GeV could be from the DD* open channel contributions. 

D. ^>'-^>(3770) mixing results 

Following the discussion of Sec. Ill C( we now investigate the energy dependence of the mixing parameter 
defined in Eq. (IT7|) . As mentioned earlier, the s-dependence of the total widths is weak for both if;' 
and V'(3770) near threshold; we thus keep the widths as constants. The g^(3?ro)DD coupling appear 
to be different in the charged and neutral channels. Here, for the purpose of investigating the energy- 
dependence of the mixing parameter, such a difference can be neglected. We then take g^uD — 9-05 
(the fitting results) and 9W377o).D£> = 12.66 (the average obtained from Ref. [281 ]) to extract the mixing 
parameter The results are presented in Fig. [5] 

The solid line denotes the mixing amplitude |£i//(s)| = |-CVv>(377o)/-CVl' while the dashed line is for 
1^-0(3770) ( s ) I = |£Vi/>(3770)/A/i(377o)l- There are clear physical meanings for these two quantities. At a 
given energy yfs, the value of |£>/>'(s)l denotes the coupling strength of the ip' component inside an initial 
-0(3770) state (propagator), while (£^(3770) denotes the strength of the -0(3770) inside an initial ip' 
state. 

The mixing parameter |£ ,y(s )l at yfs = 3.773 GeV can be related to the ip(2S)-ip(lD) state mixing 
angle via (s) | ~ | sin 0^/ | [13| for which we find 6^t ~ 5.4° . This value agrees with the one extracted in 
Ref. [j| , although it is only about half of the result estimated by Ref. [l3| . It is worth emphasizing that 
this quantity strongly depends on the coupling strength of the <ju/ d q coupling which is a consequence 
of the lineshape measurement of the cross sections near threshold. In this sense, the cross section for 
e + e~ — > DD directly provides a constraint on the mixing angle. 
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0.6- 
0.3- 



0.0- 

3.7 3.8 3.9 4.0 

Vs 

FIG. 9: The energy-dependence of the mixing parameter. The solid line is for |£<0'(s)| while the dashed line for 
1^(3770) (s)\- The vertical line labels the 0(3770) mass position. 

IV. SUMMARY AND DISCUSSION 

In this work, we investigate the reaction mechanisms for e + e~ — ¥ DD from threshold to ^fs ~ 4.3 
GeV in an effective Lagrangian approach. We find that the cross section lineshape near threshold is very 
sensitive to the presence of ip' . By fitting the cross sections from BES Q and Belle [ll[, we succeed in 
extracting resonance parameters for ip', -0(3770), A"(3900), 0(4040). The shows that the coupling gu/jjjj 
is consistent with those values given by other analyses. We show that it is important to have a reliable 
determination of this quantity, not only to understand the lineshape of the cross section at the mass of 
0(3770), but also to provide a probe for the 0'-0(377O) state mixing. 

The analysis also suggests a significant difference between 3^(3770)130 rjo and 3^(377o)_d+l>- > which is 
quite different from the PDG averaged values [28]. In order to further clarify this, a precise measurement 
of both the D°D° and D + D~ cross sections is strongly recommended, and BES-III would have great 
advantages regarding this issue (34| . 

We also study the DD* open channel effects on the e + e~ — > DD cross sections. We find that we cannot 
eliminate the fact that the enhancement at 3.9 GeV in the Belle data is due to the DD* open channel 
contributions. Further investigation of this mechanism with more accurate data would be needed. 
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